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Abstract. We classify and expose all the gradient Ricci solitons on complete surfaces, 
open or closed, with curvature bounded below, and possibly with a discrete set of cone- 
like singular points that arise naturally. We give a precise qualitative description of each 
»>f-\ metric in terms of a phase portrait, that is the most accurate description for all cases that 

y—{ do not admit an explicit expression in terms of elementary functions. Our classification 

C*^ contains examples of smooth and conic solitons that were not described in the classic 

O^ literature. We add some visual embeddings in R'^ for aesthetics. 

U 

< 

Cf) 1. Introduction 

Gradient Ricci solitons are special self-similar solutions to the Ricci flow. In [8j, R. 
Hamilton developed the Ricci flow theory for surfaces (using there a normalized version 
of Ricci flow). Hamilton proved that all closed smooth surfaces with positive curvature 
are convergent under the Ricci flow to a gradient Ricci soliton, and proved that the only 
^S^ solitons on a closed smooth surface are those of constant curvature ([8], Thm 10.1). B. 

5^ Chow subsequently was able to remove the positive curvature hypothesis in [5]. In the 

same original work of Hamilton, it is described an open gradient soliton with nonconstant 
curvature known as the steady cigar soliton. Also, in the course of the proof of Theorem 
10.1 of [8j Hamilton found some solitons on the topological sphere with cone- like singu- 
larities, and Hamilton stated that "the other solutions we have found exist on orbifolds" . 
The study of the Ricci flow converging to these orbifold solitons was carried out by L.-F. 
Wu [9], Chow and Wu [7j and Chow [3]. This brings to all 2-orbifolds a natural metric, 
\^ including those orbifolds that do not admit a constant curvature metric (footballs and 

/ teardrops). 

^^ On the other hand, a very useful fact in dimension two is that all nonconstant curvature 

C^ gradient solitons admit a nontrivial Killing vector field and have a rotational symmetry 

^^ (see [2] pp. 241-242 and[3]). This allows one to pick polar coordinates and set a single 

K>" ODE for the soliton metric. The cigar soliton appears as an explicit solution for the steady 

k> case. See [6] for further reference. 

V^ In this paper we gather and re-order these results, and we apply a phase portrait analysis 

for the ODE associated to the soliton metric. We obtain a complete classification of all 
two-dimensional gradient solitons on a surface. By "surfaces" we will mean open or closed 
topological surfaces, endowed with a complete riemannian metric, smooth everywhere 
except possibly in a discrete set of cone- like singular points, that will arise naturally. For 
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geometric considerations, on the theorem statements we will discard all riemannian metrics 
without a lower bound of the gaussian curvature, although we will find these examples in 
the course of the proofs. This solves some questions proposed in [6], p. 51, and describes 
some two dimensional solitons (smooth and conic) that didn't exist in the classic literature. 
For the sake of completeness, we will include some known results with their proofs, so this 
exposition is self-contained. The main theorem is: 

Theorem 1.1 (Main theorem). All gradient Ricci solitons on a surface, smooth every- 
where except possibly on a discrete set of cone-like singularities, complete, and with cur- 
vature bounded below fall into one of the following families: 

(1) Steady solitons: 

(a) Flat surfaces. 

(b) The smooth cigar soliton. 

(c) The cone-cigar solitons of angle a G (0, +00). 

(2) Shrinking solitons: 

(a) Spherical surfaces. 

(b) Teardrop and football solitons, on a sphere with one or two cone points. 

(c) The shrinking flat gaussian soliton on the plane. 

(d) The shrinking flat gaussian cones. 

(3) Expanding solitons: 

(a) Hyperbolic surfaces. 

(b) The a^-cone solitons, with a cone point of angle /? > and an end asymptotic 
to a cone of angle a > 0. 

(c) The smooth blunt a-cones. 

(d) The smooth cusped a-cones in the cylinder, asymptotic to a hyperbolic cusp 
in one end and asymptotic to a cone of angle a > in the other end. 

(e) The fiat-hyperbolic solitons on the plane, that are universal coverings of the 
cusped cones. 

(f) The expanding flat gaussian soliton on the plane. 

(g) The expanding flat gaussian cones. 

Each family of solitons is described in the corresponding section. In section [2] we recall 
briefly the definition of Ricci solitons in dimension 2, and their properties of symmetry in 
the case of nonconstant curvature, that lead to a first-order ODE system. In section [3] we 
enumerate the closed solitons with constant curvature. In sections |4j [5] and [6] we study the 
ODE system in the steady, shrinking and expanding cases, respectively. These three parts 
combined prove the Main Theorem. Finally, in appendix |A] we bring a gallery of solitons 
embedded into M'^, drawn with Maple. 

During the final steps of writing this paper, we found that very recently J. Bernstein and 
T. Mettler have independently posted a paper yj, that exposes a very similar classification 
for the smooth two dimensional gradient Ricci solitons, using a different method for the 
analysis of the ODEs. 

Acknowledgements. The author was partially supported by Feder/Mineco through the 
Grant MTM2009-0759. He is indebt to his advisor, Joan Porti, for all his guidance. He 
also wishes to thank Gerard Besson for the invitation to visit Grenoble in fall 2011, and 
the possibility of exposing parts of this work in the Autrans seminar in September 2011. 
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2. Gradient Ricci solitons on surfaces and Rotational symmetry 

In this section we recall the basics of Ricci sohtons and their properties of symmetry. 
See p] for an extended introduction. 

A Ricci flow is a PDE evolution equation for a riemannian metric g on a smooth 
manifold M, 

(1) |5(i) = -2Ric,(,). 

A Ricci soliton is a special type of self-similar solution of the Ricci flow, in the form 

(2) g{t) = c{t) 4>:{go) 

where for each t, c(t) is a constant and (pt is a diffeomorphism. So, g{Qi) = go, c(0) = 1 and 
(pQ = id. The family (j)t is the flow associated to a (maybe time-dependent) vector field 
X{t); and in the case when this vector field is the gradient field of a function, X = grad/, 
the soliton is said to be a gradient soliton. In this case, differenciating the soliton equation 
([2]) and evaluating at t = 0, we get 

d 
-2RlCg(o) = ^ ^_o5(0 = c(0)(/)o5'o + c(0)£grad/</'S(5o) 

= 6^0 + 2 Hessg(o) / 

where e = c(0). The soliton is said to be shrinking, steady or expanding if the constant e 
is negative, zero or positive respectively. This constant can be normalized to be —1, or 
-|-1 respectively, being this equivalent to reparameterize the time t. Therefore, a gradient 
Ricci flow is a triple {M.,g, f) satisfying 

(3) Ric + Ressf + ^g = 0. 

In the 2-dimensional case we have Ric = ^g, hence the soliton equation becomes 

(4) Hess/ + ^(i? + e)5 = 0. 

This equation only makes sense on a smooth riemannian surface. However, we will 
allow some cone-like singularities for the surface, namely points such that admit a local 
coordinate chart in the form 

dr'^ + h{r,ef d9^ 

for some smooth h : [0, 6) x M/27rZ -> M such that h{0, 9) = and ^ = ^ needs not 
to be 1. The value a is the cone angle at this point. To make this point smooth, it is 
required that a = 2tt and |^ = for all /c G N (see [6] p. 450). 

We will use some properties of two-dimensional gradient solitons to turn the tensor 
equation Q into a much simpler first order vector ODE that will allow a subsequent 
qualitative analysis. The main property we will use is that there exists a Killing vector 
field (given by a rotation of grad/) over the smooth part of M. The associated line flow 
of this field is a one-parameter group acting globally by isometries, this group must be 
S^ and so Ai is rotationally symmetric (this argument is from [2] and |3j). We define 
then a rotationally symmetric polar coordinate chart on Ai , and the analysis of the local 
expression of the soliton equation will give us the ODE system that satisfies the explicit 
metric over M. 
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Let J : TAi — ?• TAi be an almost-complex structure on M, that is, a 90° rotation on 
the positive orientation sense, so 

j2 = -Id , g{X, JX) = VX G XX. 

Lemma 1. Some basic properties of J are 

(1) g{JY,Z) = -g{Y,JZ), 

(2) J anticommutes with b: b(JX) = — J(bX), 

(3) J commutes with V. 

Proof. First property is elementary, 

= g(y + Z,J(y + Z)) = g{Y,JY)+g{Z,JY)+g{Y,JZ) + g{Z,JZ) 

= giZ,JY)+g{Y,JZ). 

A warning about the notation: if b{X) = uj is a 1-form then J : T*A4 — ;■ T*A4 is defined 
as a; I-)- Jw where Jci;(T¥) = a;(J(H^)). 
So, the second statement is just 

\){JX)(W) = g{JX, W) = -g{X, JW) = -\}{X){JW) = -J{bX){W). 

For the last statement, let X £ XAi, then {X, JX} form a basis of TAi. Let us see 
that J(VX) and V( JX) have the same projections over the basis. Since g{X, JX) = 0, 
taking covariant derivatives we get 

g{VX,JX)+g{X,V{JX)) = 

which implies 

g{J{VX), -X) + g{X, V{JX)) = 0. 

Again, differenciating g{X, X) = g{JX, JX) we obtain 

2g{X,VX) = 2g{V{JX),JX), 

which implies 

g{JX,J{VX))=g{JX,V{JX)). 

D 

This construction shows that a two-dimensional gradient soliton admits a Killing vector 
field. 

Lemma 2 ([2] p. 241). The vector field J(grad/) is a Killing vector field. 

Proof. Recall that a Killing vector field W is such that its line flow is by isometries, or 
equivalently, the metric tensor g is invariant under the line flow, that can be expressed in 
terms of the Lie derivative as Cwg = 0. Recall also that 

Cwg {y, z) = VL^(y, z) + Vw(z, y) 
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where uj = \>{W) = g{W, •). Then 

^J(grad/)9(i",^) = V(b(J(grad/)))(y,Z)+V(b(J(grad/)))(Z,y) 
= VJbgrad/(y,Z) + VJbgrad/(Z,y) 

= vjv/(y,z) + vjv/(z,y) 
= jvv/(y,z) + jvv/(z,y) 

= V^f{JY,Z) + V^f{JZ,Y) 

= ^iR + e)g{JY,Z) + ^{R + e)g{JZ,Y) 



~{R + e)(^g{JY,Z) + g{JZ,Y)) =0. 



Note that Vf = df = b(grad/), and again J(V/) is the 1-form A H> Vf{JA) and JiV^f) 
is the 2-covariant tensor field {A, B) ^ V^fljA, B). U 

The Kilhng vector field may be null if the grad / field itself is null, otherwise, the surface 
admits a symmetry. 

Lemma 3. Let (Ai,g,f) be a gradient Ricci soliton on a surface. Then, at least one of 
the following holds: 

(1) M. has constant curvature. 

(2) M. is rotationally symmetric (i.e. admits a S^ -action by isometrics). 

(3) A4 admits a quotient that is rotationally symmetric. 

Besides, if the surface has not constant curvature, no more than two cone points may exist. 

Proof. We adapt the argument for the smooth closed case from [3j. We will discuss in 
terms of grad/. If grad/ = 0, then by the soliton equation Q we have R = —e and the 
curvature is constant. Let us assume then that / is not constant everywhere. Therefore 
J(grad/) is a nontrivial Killing vector field and its line flow, (j)t, is a one-parameter group 
acting over A4 by isometries. 

Suppose that grad / has at least one zero in a point O £ M.. This is the case of closed 
smooth surfaces. The point O is a zero of the vector field J(grad/), so it is a fixed point 
of (pt for every t. Then, (j)t induces (p^ acting on TqM by isometries of the tangent plane, 
so we conclude that the group {(pt} is S^ acting by rotations on the tangent plane. Via the 
exponential map on O, the action is global on Ai and therefore the surface is rotationally 
symmetric. 

Suppose now that grad / has no zeroes but the surface contains a cone point P. Then 
the flowlines of (pt cannot pass through P, because there is no local isometry between a 
cone point and a smooth one. So this point P is fixed by (pt for every t and, via the 
exponential map, (pt induces (pt acting on CpM the tangent cone (space of directions) on 
P. Again, a continuous one-parameter subgroup of the metric cone CpAi must be the S^ 
group acting by rotations. Besides, if other cone points were to exist, these should also 
be fixed by the already given §^ action. This implies that no more than two cone points 
can exist on M, for otherwise the minimal geodesies joining P with two or more conical 
points would be both fixed and exchanged by some S^ group element. Note that in the 
case of two cone points, these need not to have equal cone angles. 
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Finally suppose that grad / has no zeroes and the surface has no cone points. Then the 
surface is smooth and the flowlines of grad / are all of them isomorphic to M (no closed 
orbits can appear for the gradient of a function) and foliate the surface. The action of (j)t 
exchanges the fibres of this foliation. The parameter of </>( is t G S^ or t G M. In the first 
case, S^ is acting on M and it is rotationally symmetric. In the second case, A^ = M^, 
and the flowline (pt of the Killing vector field induces a Z-action by isometrics by 

X I— ;■ 4>i{x) 

that acts freely on Ai since no point is fixed by (j)t for any t 7^ (if (t)t{p) = P, then 
all fibres are fixed and every point in each fibre also is, so (j)t = id)- Then the quotient 
by this action is topologically M./ ^ = M x S^ and is rotationally symmetric. We will 
find nontrivial examples of this solitons as cusped expanding solitons and their universal 
coverings. D 

The fact of being rotationally symmetric allows us to endow M with polar coordinates 
(r, 9) such that the metric is given by 

where r G / C M is the radial coordinate, and 6 G M/27rZ is the periodic angular co- 
ordinate. The function h{r) does not depend on 9 because of the rotational symmetry; 
and similarly, the potential function only depends on the r coordinate, since grad/ is a 
radial vector field. Surfaces not rotationally symmetric but with a rotationally symmetric 
quotient also admit these coordinates, changing only G M. 

Lemma 4. Given the polar coordinates (r, 0) G M x R/27rZ and the metric in the form 
g = dr"^ + /i^(r) dO"^ , the gaussian curvature (which equals half the scalar curvature) is 
given by 

'< = - = =f. 

2 h ' 

and the hessian of a radial function f{r) is given by 

Hess/ = /" dr^ + hh'f d9^. 

Proof. It is a standard computation. Covariant derivatives are given by 

h' 
Vg^dr = Vg^de = j-de ^dgde = -hh'dr 

Then we contract twice the curvature tensor R{X,Y)Z = Vx(VyZ) — Vy(Vx-^) — 
V[x,Y]^ ^^^ ^^^ scalar curvature, and apply Hess/(X, Y) = X{Y{f)) — {VxY){f) for the 
hessian. □ 

On that rotationally symmetric setting, the soliton equation becomes 

Hess/ + ^(i? + 6)5=(r-^ + 0^r2+(Mr+(-y + 0/i^)d0^ = O, 
which is equivalent to the second order ODEs system 



(5) 



r-^ + i = 
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We combine both equations to obtain 

f h' 
and integrating this equation, 

Inf =lnh + C 
so 

(6) /' = ah 

for some a > 0. Hence, substituting on the system we obtain a single ODE, 

(7) h" - ahh' --h = 0. 

We summarize the computations in the following lemma. 

Lemma 5. Let {Ai,g,f) be a gradient Ricci soliton on a surface with nonconstant cur- 
vature. Then Ai admits coordinates (r, 9), with r € / C R and 6 ^ S^ or 9 gR, such that 
the metric takes the form g = dr"^ + h'^{r) dO"^ for some function h = h{r) satisfying (JTl), 
and the potential is f = f{r) satisfying ^. 

Setting h' = u, the second order ODE (JTJ) becomes a vector first order ODE 

h' = u 



^^^ ^ u' = {au+^)h. 



The solutions to system ([8]) (and equation M) are functions h{r) that define rotationally 
symmetric metrics on the cylinder {r,9) G M x S^. This cylinder may be pinched in one 
or both ends, thus changing the topology of the surface. The pinching appears as zeros of 
h. Closedness condition of the surface is equivalent to the boundary conditions 

h{0) = and h{A) = 

for some ^4 > such that h{A) = 0. In this case, one or two cone angles may appear, 

h'{0) = ^ and h'{A) = -2 

where ai and 02 are the cone angles. Smoothness conditions would be h'{0) = 1 and 
h'{A) = — 1. We shall study the system ([s]) for steady, shrinking and expanding solitons 
to obtain a complete enumeration of gradient Ricci solitons on surfaces of nonconstant 
curvature. 

Incidentally, it is interesting to note some geometric interpretations of the functions / 
and h. 

Lemma 6. Let {Ai,g,f) be a rotationally symmetric gradient Ricci soliton in dimension 
2, then 

grad/ = f'{r)dr = ah{r)d.r, 

fir) 
J(grad/) = -ri-^dg = ade, 
n{r) 

K = -'^ = -(ah'+'- 
h V 2 
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If p G M is a (smooth or conic) center of rotation. Then, using the distance to p as the 
r- coordinate, 

h[r) = — Perinieter(Disc(r)), 

27r 

/(r) = ^Area(Disc(r)) + /(0), 

where Disc(r) is the disc with radius r centered at p. On the other hand, if there is no 
center of rotation, then 

Hro) = 7^Length(L), 

ZTT 

/(n))-/(ro) = |^Area(S), 



where L is the level set {r = r^} and B is the annulus hounded by the two level sets 
{r = ro} and {r = ri}. 

3. Closed solitons of constant curvature 

Before looking for the specific nontrivial steady, shrinking and expandig sohtons, in this 
section we rule out the constant curvature solitons that also are rotationally symmetric. 
Following [3], if we look for rotationally symmetric closed smooth solitons (/i(0) = h{A) = 
and h'{0) = —h'{A) = 1), we can show that there is no other function / than a constant 
one. For we multiply equation Q by h' to get 

hh' 

h'h" - ah{h'f + — = 



and integrate on [0, A\ to obtain 




{h'f 
2 


- a / h{h'fdr + — 
-^o ^ 


A 



and since /i(0) = h{A) = 0, and /i'(0) = -h'{A), 




f-A 

= -a / h{h'fdr 
Jo 





which is impossible unless a = 0. This is indeed the case when /' = 0, there is no gradient 
vector field, no Killing vector field, constant curvature and the soliton is a homothetic fixed 
metric. Note that if there is no vector field, there is no need to be rotationally symmetric, 
thus one can have constant curvature surfaces of any genus. Therefore we have seen the 
following lemma. 

Lemma 7. The only solitons over a compact smooth surface are those of constant curva- 
ture. 
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More generally, rotationally symmetric closed solitons with two equal angles satisfy 
h'{0) = —h'{A) = ^ and the same argument applies. In this case, equation ([T]) turns into 

h" - |/i = 
that can be explicitly solved. For e = 1 the solution is 

h{r) = cie"/^ + C2e-'^/^ 

but the closedness condition /i(0) = h{A) = implies ci = C2 = 0. Thus there are no 
expanding solitons with two equal cone points besides the constant curvature ones. 

For e = 0, the solution is h{r) = cir+C2, that can't have two zeroes unless h = 0. Finally, 
for e = — 1 the solution is h{r) = ci sin(r/\/2) + C2 cos(r/\/2), and by the closedness C2 = 0. 
This metric is locally the round sphere. We have therefore seen the following. 

Lemma 8. The only solitons over a compact surface with two equal cone points are shrink- 
ing spherical surfaces. 

Up to now, we have examinated all possible cases with a = (equivalently, with / 
constant and with constant curvature). Thus, we will assume henceforth that a ^ and 
/ is not constant. 

4. Steady solitons 

In this section we study the steady case (e = 0) of rotationally symmetric solitons. The 
equation (JTJ) reduces to 

(9) h" - ahh' = 

and the system ([s]) to 

(>») { t' : luH 

The phase portrait of ( 10 ) is shown in Figure [l] 

This phase portrait has a line of fixed points at {u = 0}, that account for the trivial 
steady solitons consisting on a flat cylinder of any fixed diameter (or their universal cover- 
ing, the flat plane). No other critical points are present. Only the right half-plane {h > 0} 
is needed, since we can take /i > in the metric definition. 

Every integral curve of the system lies on a parabola. This follows from manipulating 



system (10) 



ahh' = a\ — 
2 



and hence 

/l2 



u = a— + C. 



In another terminology, the function 



/i2 
H{h, u) = a u 
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Figure 1. Phase portrait of the system (10) with a = 1. 



is a first integral of the system (10). Furthermore, we can finish the integration of the 
equation 



by writting 

The solution to this ODE is 

(11) h{r) 

if C7 > 0; 

(12) h{r) = 

if C < 0; and 

(13) 



h' = a— + C 
2 



h' 



c + {VlhY 




-2C 



h{r) 



tan 



tanh 




r + D 



-aC 



r + D 



D-^r 



if C = 0. 

Now, let us examinate each type of solution. If C > 0, the parabola lies completely on 



the upper half-plane {u > 0}. The equation ( 11 ) implies that /i — t- oo for some finite value 
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of r, and hence the metric is not complete. Furthermore, the gaussian curvature of the 
metric K satisfies 

u = -Kh 

and since u is increasing on these solutions, the curvature is not bounded below. The case 
C = 1 is sometimes called the exploding soliton in the literature ([6]). 

If we look at C = 0, the parabola touches the origin of coordinates, and its right 
hand branch defines a metric on the cylinder. The value of D = j^j^ can be set so that 

D = h{0) = 1 just reparameterizing r. With this parameterization, r € (— oo, -. For r < 0, 
the function h is well defined and determines a negatively curved metric that approaches 
a cusp as r — )• — oo. However, for r G [0, -) the metric is not complete and its curvature 
tends to — oo as t — )■ -. 

a 

We look now at the case C < 0, first for the solutions lying in the lower half-plane 
{u < 0}. We can assume h{0) = 0, C = u(0), D = and r < (this means that — r is the 
arc-parameter of the meridians). All these arcs of parabolas join a point on the {h = 0} 
axis with a point on the {h' = u = 0} axis. This means that the cylinder is pinched in 
one end, and approaches a constant diameter cylinder on the other end. The metrics are 



complete on the cylindrical end, because from equation (12) /i — t- est as r — t- — oo. The 
curvature on these metrics is bounded and positive, since u and u' < are bounded. 

Some of these metrics are smooth, the particular cases of C = n(0) = /i'(0) = —1. Note 
that derivating the equation h" = ahh' and evaluating at r = one sees that all even-order 
derivatives vanish and the surface is truly C°° at this point. These are the so called cigar 
solitons. There are actually infinitely many of them, adjusting the value of a and changing 
the diameter of the asymptotic cylinder, although all of them are homothetic and hence it 
is said to exist the cigar soliton. All the other metrics have a cone point at r = 0, whose 
angle is — 27r/i'(0). 

The only remaining case to inspect is the solutions with C < lying on the upper half 
plane {u > 0}. These unbounded arcs of parabolas rise from the axis {u = 0}. We can 
assume (changing D and reparameterizing r) that r G [0, -|-oo). The metric is complete in 



r — )• +00 because of equation (12), however, these metrics fail to be complete on r = 0, 
having a metric completion with boundary S^. The curvature on these metrics is negative 
and not bounded below. 

We summarize the discussion in the following theorem. 

Theorem 4.1. The only complete steady gradient Ricci solitons on a surface with curva- 
ture bounded below are: 

(1) Flat surfaces (possibly with cone points). 

(2) The smooth cigar soliton (up to homothety). 

(3) The cone-cigar solitons of angle a G (0, -|-oo) (up to homothety) . 

Remark. There exist other steady gradient solitons, with curvature not bounded below, 
as described above. 

Pictures of a cigar soliton and a cone-cigar soliton are shown in Figures [6] and [7| 

5. Shrinking solitons 

In this section we study the shrinking solitons (e = —1), besides the round sphere and 
the spherical footballs with two equal cone angles found in section |3j When e = — 1, the 
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metric of Ai is determined by a real- valued function /i(r) satisfying the second order ODE 

h 



(14) 

or equivalently the system 

(15) 



h" — ahh' + 



0, 



h' 



u 

{au — h)h. 



The phase portrait of this ODE system with a = 1 is shown in Figure [2] This phase 
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Figure 2. Phase portrait of the system (15) with a = 1. 

portrait has a critical point at (/i, u) = (0, 0) of type center, and a horizontal isocline 
(points such that u' = 0) at the line ^ = 2a' -^^-ch curve on this /iit-plane corresponds 
to a solution h, and the intersection with the vertical axis {h = 0} are at u(0) and u{A), 
which stand for the cone angles. Indeed, only half of each curve is enough to define the 
soliton, the one lying in the {h > 0} half-plane, since we can choose the sign of h because 
only /i'^ is used to define the metric. 

All curves in the phase portrait represent rotationally symmetric soliton metrics over, 
a priori, a topological cylinder. Closed curves (that intersect twice the axis {h = 0}) 
are actually metrics over a doubly pinched cylinder, thus a topological sphere with two 
cone points, giving the so called football solitons. Open curves only intersect once the 
{h = 0} axis, and hence are metrics over a topological plane. If the intersection of any 
curve with the {h = 0} axis occurs at n = ±1, then the metric extends smoothly to this 
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point (truly C°° since derivating (14) all even-order derivatives vanish at this point). For 
instance, in Figure [2] there is only one curve associated to a teardrop soliton, namely the 
one intersecting the vertical axis at some value n(0) G (0, 2) and at u{A) = —1. There 
is also a smooth soliton metric on M^, namely the one associated with the curve passing 
through {h,u) = (0,1), and all other curves represent solitons over cone surfaces. The 
separatrix line, u = j^, represents the solution h{r) = 5^ + co, which stands for the metric 
dr^ + j^r'^dO'^ . This is a flat metric on the cone of angle ^, a cone version of the shrinking 
gaussian soliton, and we call it a shrinking gaussian cone soliton (the smooth shrinking 
gaussian soliton is the case a = |). 



Let us focus on the compact shrinking solitons. 

Lemma 9. For every pair of values < ai < a2 < 00, there exist a unique value a > 
such that the equation ([7]) has one solution satisfying the boundary conditions h'{0) = ^ 
andh'{A) = -^. 

Equivalently, the lemma asserts that there exists a value a such that the phase portrait of 
the system (pi) has one solution curve that intersect the vertical axis {h = 0} at n(0) = ^ 
andn(A) = -g. 

Proof. We can normalize the system by 

V = ah 
w = au 

so that on this coordinates the system becomes 

(16) I "; = r 1, 

This would be the same system as (15) with a = 1, which is indeed shown on Figure [2] 
The system (16) has the following first integral, 

H{v,w) = v — 2w — lii\2w — 1\ 

that is, the solution curves of the system are the level sets of H. Indeed, derivating 
H{v{r),w{r)) with respect to r, 

-p—ntVjW) = 2vv — 2w 

or 2vu — 1 

= 2vw-2 [w- -] v{l + ) =0. 

The cone angle conditions are ai = ^"'^ , 02 = — '^^ , while t'(O) = v{A) = 0. 
Thus, the function w evaluated at and A satisfies 

H{<d, w) = 2w + hi\2w - 1\ = C 

for some C &M.. This is equivalent, via 2w — 1 = —y and e*^ = k, to the equation 

(17) \y\ = ke^-^ 

(cf. [8]). Although not expressable in terms of elementary functions, this equation has 
three solutions for y, one for negative y and two for positive y (see Figure pi). The two 



positive solutions of (17) are the intersection of the exponential function e^ ^ with the 
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line ^y with slope ^. These two positive solutions are associated to a compact connected 
component of H{v,w) = C, whereas the negative solution is associated to a noncompact 
component of H that represent noncompact soliton surfaces. The two positive solutions 




Figure 3. The graphs of the exponential e^ ^ and ^|y| 



of (17) exist only when k G (0, 1) and actually these two solutions are equal when k = 1 
and the line is tangent to the exponential function at y = 1. These two solutions yi, y2 of 



(17) are therefore located on (0, 1) and (1, +oo) respectively, and can be expressed as 

2/i = l-p , 2/2 = 1 + 9 
with p,q > 0. The two cone angles, having assumed ai < 02, are then expressed as 

27rt(;(0) 27r 1 — 2/1 vr 
a a 2 a~ 

2ttw{A) 27r 1 - y2 



«! 



02 



= 27r/i'(0) = 27ru(0) : 
2TTh'{A) = -27ru{A) 



-P 



IT 



and their quotient is 



02 



p 



Let * : (0, 1) 



be the mapping 



k H> ^{k) 



p 



The function ^ is injective and the quotient ^(k) ranges from to 1 when varying k E 
(0, 1). This is proven in |8j. Lemma 10.7, we can visualize its graph in FigureHj Therefore 
for any pair of chosen angles ai < 02 there exist k = ^ 



(17) has two positive solutions yi. 



^'^), such that the equation 
2/2- This yields two values p = 1 — yi, q = 1/2 — 1, and 



finally we recover 

ai 02 

vrp nq 

This value makes the system ([8| and the equation (JTl) to have the required solutions. D 
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0.4 0.6 

k 



0.8 



Figure 4. The function ^!{k). 

We summarize our discussion in the following theorem. 

Theorem 5.1. The only complete shrinking gradient Ricci solitons on a surface 7W with 
curvature hounded below are: 

• Spherical surfaces (including the round sphere and the round football solitons of 
constant curvature and two equal cone angles). 

• The football and teardrop solitons, if M. is compact and has one or two different 
cone points. These cone angles can be any real positive value. 

• The smooth shrinking flat gaussian soliton. 

• The shrinking fiat gaussian cones. 



6. Expanding solitons 



We end our classification with the expanding solitons (e = 1). The equation ph and 
the system ([8|) are in this case 



(18) 
and 
(19) 



h — ahh -\ — = 
2 



h' 



[au 



+ l)h. 



The phase portrait of (19) is shown in Figure pi We can rescale the system (19) with the 
change 

V = ah 

w = au 

so that on this coordinates the system becomes 



(20) 



V 

w' 



w 



[W + 2Jf 
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Figure 5. Phase portrait of the system (19) with a = 1 



which is exactly the system (19) with a = 1. The phase portrait of this system is indeed 



shown on Figure [5| We will study the trajectories of the normalized system and next we 
will discuss the geometrical interpretation of each trajectory. 

The system ( |20[ ) has a critical point [v' = w' = 0) at (0, 0). It has an horizontal isocline 
[w' = 0) at the line L = {w = — gli that is also an orbit solution, and hence no other 
trajectory can cross it. The vertical axis {v = 0} is also an horizontal isocline. The 
horizontal axis {w = 0} is, on the other hand, a vertical isocline {v' = 0). 



The linearization of the system (20) at the critical point (0,0) is 



w 





w + 



V 

w 



The matrix of the linearized system at the critical point is 



, that has determinant 



— 2 < and hence the critical point is a saddle point. The eigenvalues of this matrix are 
4= and — ^, with eigenvectors respectively 



72 
1 



and 



-V2 
1 



These eigenvectors determine the two principal directions of the saddle point, from which 
four separatrix curves are emanating. 



The system (20) has the following first integral, 

,2 



H{v, w) 



V 



2w + \n\2w + l\ 
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that is, the solution curves of the system are the level sets of H. Indeed, derivating 
H{v{r),w{r)) with respect to r, 

-H[v,w) = 2vv — 2w + 



9r ' 2w + 1 

= 2vw-2 (w + -]v(l ) =0. 

V 2j \ 2w + l) 

From the system, and more apparently from the first integral, it is clear that the phase 
portrait is symmetric with respect to the axis {v = 0}. We will only study then the 
trajectories on the right-hand half-plane {u > 0}. Actually this restriction agrees with the 
geometric assumption of /i > 0. 

Let us consider a trajectory passing through a point in the quadrant {v > 0,w > 0}. 
Then v' > and w' > and hence the curve moves upwards and rightwards. More 
carefully, as soon as f > 5i > and w > 62 > 0, both derivatives are bounded below away 
from zero, v' > 5i > and w' > 62 > 0, and therefore v and w tend to -|-oo. We can 
further evaluate the asymptotic behaviour from the first integral, 

v"^ ln|2tf -Ml _ C 

2w 2w 2w 

We take the limit as r — )■ -|-oo and since w — t- -|-oo, we get that 

lim ^ = 1, 
r^+oo 2w{r) 

so the orbit approaches a parabola (the same parabolas of the steady case and asymptotic 
on the shrinking case). 

We now inspect the separatrix S emanating (actually sinking) from the critical point 
at the direction (v2, ~1)- The associated eigenvalue is — v2 and hence the trajectory is 
approaching the saddle point (hence the sinking). The curve S lies in the w' > region, 
and cannot cross the horizontal isocline L. Therefore, the separatrix when seen backwards 
in r must be decreasing and bounded, and hence must approach a horizontal asymptote. 
This asymptote must be L, since if the trajectory were lying in the region w' > 6 > 
for infinite time, it would come from w = —00, which is absurd since it cannot cross the 
isocline L. Therefore, over this separatrix S, v ^ -|-oo and w ^ —^ as r — )■ —00. 

Any trajectory lying over S will eventually enter in the upper right quadrant, and 
hence approach asymptotically the previously mentioned parabolas. This is clear since w' 
is positive and bounded below away from zero. 

Let us study the trajectories below S. All these curves intersect the axis {u = 0}, and 
we can consider the origin of the r coordinate as such that the intersection point with the 
axis occurs at r = 0. Then, the region of the curves parameterized by r < lies in the 
V > half-plane. Since the curves are below S, they lie in the lower right quadrant and 
hence v' < and f— J'-l-ooasr— t-— 00. If the curve lies over L, then w' > 0, and if it lies 
below L, then w' < 0. This means that the isocline is repulsive forward in r and attractive 
backwards in r. Therefore any curve lying below S will have an asymptote as r — t- —00 
and as before this must be the isocline L, that is, v — )■ -|-oo and w ^ —^ as r — )■ —00. 

Let us comment about the domain of r. We have stated that the trajectories below 
S are parameterized for r G (— cx),0], although a priori it could be r G (— M, 0] for some 
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maximal M (and hence v — t- +oo as r — t- — M, and these would represent noncomplete 
metrics). This is not the case, since the trajectories are approaching v' = — j, and hence 
v' is bounded (|f'| < 1 for r less than some ro < 0), so f cannot grow to +oo for finite 
r-time. 

Finally, let us study the separatrix S itself. This curve is parameterized by r G M, and 
{v,w) — )• (0,0) as r — )• +oo and {v,w) — )• (+oo,— ^) as r — )• — oo (this follows from the 
Grobman-Hartman theorem in the end near the saddle point, and from the asymptotic 
L on the other end). We can give a more detailed description of the asymptotics. As 
r — )• — oo, we know that w — >■ — ^, this is 

lim — J- = 1. 

r— >— oo — — 
2 



Then, applying the I'Hopital rule, 



lim — T— = 1, 

r— >~oo — — T 
2' 



or v{r) ~ — 2^ as r — )■ — oo. This is valid for all the trajectories asymptotic to the 
horizontal isocline. Similarly, as r — t- +oo, we know that f,i(; — t- 0, but furthermore we 
know that their quotient tends to the slope of the eigenvector determining the separatrix, 
i.e. 

lim — = lim — = —^, 

r-5>+oo w r-s>+oo v -y/2 



which is to say 

Then, by I'Hopital rule. 



lim (Inv)' = lim — = — v2. 

•— >+oo r— >+oo V 



,. (Inw) Int; 

lim ^ = lim -j=- = 1. 

'-s>+oo — y2 r->+oo — -y/2r 



this is, v{r) ~ e"^^'' as r — t- +oo. 

At this point we have established all the important features of the phase portrait in 
Figure^ With the unnormalyzed system, in coordinates (/i, n), the phase portrait is just 
a scaling of the one in Figure [5] by the factor a, and thus the horizontal isocline L is at 
{u = —2^}. We now give the geometric interpretation of each trajectory. 

All the trajectories above S have solutions with unbounded positive w. This means 
that the gaussian curvature of the associated metric K = —{au + \) = —{'w + \) is not 
bounded below, and we will discard them. 

All the trajectories below the separatrix S have bounded w and therefore bounded 
curvature on the associated metric. More specifically, the curves above the isocline L will 
give metrics with negative curvature, and curves below L will give metrics with positive 
curvature. These curves will intersect the {h = 0} axis at 6 < 0, and the associated metric 
will have a cone point of angle 

/3 = -2-Kh 

at the point of coordinate r = 0. On the other end, the function h{r) is asymptotic to 
— 2^r (recall that the parameter is r G (— oo,0]) and the metric will be asymptotic to the 
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wide part of a flat cone of angle 



IT 

a = —. 
a 



We call these solitons the aP-cone solitons. These solitons have positive curvature if a < /3 
(6 < ^) and negative if a > f3 {b > ^). In the case a = (3 [h = ^) we are in the case 
of the isocline L. This line h' = u = —i^^, has as solution the parameterization 

h{r) = ^^r + C 

which represents a flat expanding gaussian cone soliton, with cone angle -. The special 
case a = 2 yields a smooth metric at r = 0, thus we have a flat metric on the plane known 
as the flat expanding gaussian soliton. 

Other remarkable cases are those with (3 = 2-k (b = —1), because the cone point at the 



apex is now blunted and the surface is smooth (we can check from equation ( 18 ) that all 
even-order derivatives vanish at r = 0), we call them the blunted a-cone solitons. The 
angle a may be less or greater than 2tt (but in the later case it cannot be embedded 
symmetrically in M^). 

We interpret now the separatrix S. This is the limiting case as the angle (3 tends to zero. 
In this case the parameter r is not on (— oo, 0] but on the whole M and thus h{r) defines a 
smooth complete metric on the cylinder. As r — )• +c«, the function h{r) is asymptotic to 
_g-v2r^ that defines a hyperbolic metric of constant curvature —2. This hyperbolic metric 
on a cylinder is called a hyperbolic cusp. The separatrix L represents a soliton metric that 
approaches the thin part of a hyperbolic cusp in one end, and the wide part of a flat cone 
on the other. There is still freedom to set the angle a, and we call these the cusped a-cone 
solitons. 

Finally, there is still one more family of two dimensional gradient solitons, namely the 
universal cover of the cusped a-cones. These solitons are metrics on M? locally isometric 
to the cusped cones. These solitons are not rotationally symmetric, but translationally 
symmetric, i.e. there is not a S^ group but a M group acting by isometries. The plane M^ 
with any of this metrics has a fixed direction (given by grad /) such that a straight line 
following this direction (that is also a geodesic of the soliton metric) transits gradually 
from a region of hyperbolic curvature on one end to a region of flat curvature on the 
other. Any translation on the direction perpendicular to grad / (this is, in the direction of 
J(grad(/) ) is an isometry on these metrics. We call these flat-hyperbolic soliton planes. 

We summarize our discussion in the following theorem. 

Theorem 6.1. The only complete expanding gradient Ricci solitons on a surface Ai with 
curvature bounded below are: 

• Hyperbolic surfaces (possibly with cone angles). 

• The af3-cone solitons, for every pair of cone angles a,/3 > 0. 

• The smooth complete blunt a-cones, that are aj3-cones with (3 = 2tt. 

• The smooth expanding flat gaussian soliton. 

• The expanding flat gaussian cones. 

• The smooth cusped a-cone in the cylinder. 

• The flat-hyperbolic solitons on the plane, that are universal coverings of the cusped 
cones. 
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Appendix A. Gallery of embedded solitons 

Just for aesthetics, we can embed some of the solitons we described into M'^ and visuahze 
them numerically as surfaces with the inherited metric from the ambient euclidean space. 
If we want to keep the rotational symmetry apparent, however, we cannot embed into 
M^ a cone point of angle greater than 27r, and we can't embed a rotational surface whose 
parallels have length L = 2tt h{R) if R < i^. 

In order to do this, we use the metric in polar coordinates (r, 9) G [0, A] x [0, 27r], 

dr^ + h{rf d9^. 

We recall that r is the arc parameter of the {0 = est} curves (meridians), and that 
the {r = est} curves (parallels) are circles of radius h{r) parameterized by 6* G [0,27r]. 
Therefore, we can use the h, 6 as polar coordinates on the plane, and find an appropriate 
third coordinate z (height). When we put the stacked parallels of radius h{r) at height 
z{r), we obtain a rotational surface whose meridians have length parameter r. Thus, 



dr^ = dh^ + dz^ 



or equivalently 

which defines z = z{r) as satisfying 



dz2 _ dh'^ 
dr^ dr"^ 



{z'f = 1 



.2 



U 



with the convention that h' = u. Hence, to obtain an embedded surface satisfying the 
soliton system ([8]) it is sufficient to integrate the first order vector ODE 

h' = u 

u' = {au + |)/i 

z' = ^fY~^ 

with initial conditions /i(0) = 0, u(0) = 6, 2;(0) = 0. Once obtained a numerical solution 
for /i(r), u{r^ and z{r\ we can fix a value ^ > and then plot the set of points 

{(/i(r)cos6',/i(r)sin6',z(r)) G M''^ | r G [0, yl] , 6^ G [0,27r]}. 

Next we show some embedded solitons. These were obtained with the following Maple 
code: 

> epsilon:=l; a:=l; b:=-l; A:=10; 

> sys:= diff (h(r) ,r)=u(r) , diff (u(r) ,r)=(a*u(r)+epsilon/2)*h(r) , 
diff (z(r) ,r)=sqrt(l-u(r)~2) : 

fns:= dsolveC {sys , h(0)=0, u(0)=b, z(0)=0}, numeric, 

output=listprocedure) : 

hh:=rhs(fns[2]) : uu:=rhs(fns [3] ) : zz:=rhs(fns [4] ) : 

> plot3d( [hh(r)*cos(theta) , hh(r)*sin(theta) ,zz(r)] , 

r=0 . . A , theta=0 . . 2*Pi , scaling=constrained , grid= [40 , 40] ) ; 
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Figure 6. A cigar soliton (e = 0, a = 1, 6 = —1). 



Figure 7. A cone-cigar soliton with cone angle 180° (e = 0, a = 1, 6 = —0.5). 




Figure 8. A football soliton with cone angles 108° and 183.38° (e = -1, 
a = 1, 6 = 0.3, yl = 4.56). 




Figure 9. A teardrop soliton with cone angle 169.36° (e = —1, a = 0.8, 
b= -1, ^ = 4.68). 
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Figure 10. An a/3-cone soliton with cone angles a = 240° and (3 = 90° 
(e = 1, a = 0.75, b = —0.25). Note that the curvature is negative since 

a > /3. 




Figure 11. An Q/3-cone sohton with cone angles a = 180° and /3 = 306° 
(e = 1, a = 1, 6 = —0.85). Note that the curvature is positive since a < j3. 




Figure 12. A blunt a-cone soliton with asymptotic cone angle a 
(e = 1, a = 1, 6 = -1). 



180° 




Figure 13. A cusped a-cone soliton with asymptotic cone angle a 
(e = 1, a = 1, separatrix S {h ^ 0)). 



180° 
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